First forbidden /3-decay rates for neutron-rich nickel isotopes in stellar matter 
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In astrophysical environments, allowed Gamow- Teller (GT) transitions are important, particularly 
for /3-decay rates in presupernova evolution of massive stars, since they contribute to the fine-tuning 
of the lepton-to-baryon content of the stellar matter prior to and during the collapse of a heavy star. 
In environments where GT transitions are unfavored, first-forbidden transitions become important 
especially in medium heavy and heavy nuclei. Particularly in case of neutron-rich nuclei, first- 
forbidden transitions are favored primarily due to the phase-space amplification for these transitions. 
In this work the total /3-decay half-lives and the first-forbidden(FF) /3-decay rates for a number of 
neutron-rich nickel isotopes, 72_78 Ni, are calculated using the proton-neutron quasi-particle random 
phase approximation (pn-QRPA) theory for the first time in stellar matter. For the calculation of 
the /3-decay half-lives both allowed and first-forbidden transitions were considered. Comparison of 
the total half-lives is made with measurements and other theoretical calculation. The pn-QRPA 
results agree reasonably well with experiments. 

PACS numbers: 21.60.Jz, 23.40. Bw, 23.40.-s, 26.30. Jk, 26.50.+X 



I. INTRODUCTION 



Reliable and precise knowledge of the /3-decay for neutron-rich nuclei is crucial to an understanding of the r-process. 
Both the element distribution on the r-path, and the resulting final distribution of stable elements are highly sensitive 
to the /3-decay properties of the neutron-rich nuclei involved in the process Q, Q] ■ There are about 6000 nuclei between 
the j3 stability line and the neutron drip line. Most of these nuclei cannot be produced in terrestrial laboratories and 
one has to rely on theoretical extrapolations in respect of beta decay properties. In neutron-rich environments electron 
neutrino captures could not only amplify the effect of /3-decays but the subsequent i/-induced neutron spallation can 
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also contribute towards changing the r-abundance distribution pattern 3]. Correspondingly reliable predictions of 
/3-decay for neutron-rich nuclei are considered to be very important for r-process nucleosynthesis. 

The weak interaction rates are the important ingredients playing a crucial role in practically all stellar processes: 
the hydrostatic burning of massive stars, presupernova evolution of massive stars,and nucleosynthesis (s-, p-, r-, rp-) 
processes (see, for example, the seminal paper by Burbidge and collaborators [J]). For densities p < 10 n g/cm 3 , 
stellar weak interaction processes are dominated by Gamow- Teller (GT) and, if applicable, by Fermi transitions. For 
nuclei lying in the vicinity of the line of stability, forbidden transitions contribute sizably for p > lO^g/cm 3 when 
the electron chemical potential reaches values of the order of 30 MeV or more 5]. 

During presupernova stage of stellar evolution, electron capture dominates the weak interaction processes. However, 
after silicon depletion in the core and during the silicon shell burning, /3-decay competes temporarily with electron 
capture and further cools the star. Unlike electron capture, the /?-decay at silicon burning stage increases and this 
bears consequences because /3-decays are additional neutrino source and add to the cooling of the stellar core and a 
reduction in stellar entropy. This cooling of stellar core can be quite efficient as often the average neutrino energy in 
the involved /3-decays is larger than for the competing electron capture processes. Consequently, the /3-decay lowers 
the stellar core temperature after silicon shell burning. For a discussion on the fine-tuning of lepton-to-baryon ratio 
of stellar matter during presupernova evolution see Ref. 6| . As the density of the stellar core increases the allowed 
/3-decay becomes unimportant and hindered due to the fact that the increasing Fermi energy of electrons blocks the 
available phase space for the electron to be produced in the /3-decay. (For a detailed discussion of the role of weak 
interaction in the presupernova evolution of massive stars see Ref. p}.) In case of neutron-rich nuclei, first-forbidden 
/3-decay may become important due to the enlarge phase space for these transitions 8|, [9]. 

Reliable, quantitative estimates of /3-decay half-lives of neutron-rich nuclei are needed in astrophysics for the 
understanding of supernova explosions, and the processes of nucleosynthesis, particularly the r-process. The /3-decay 
half-lives are also needed for the experimental exploration of the nuclear landscape at existing and future radioactive 
ion-beam facilities. The calculation of /3-decay half-lives in agreement with experimental results has been a challenging 
problem for nuclear theorists. 

Because of the scarcity of experimental data, majority of the /3-decay rates of the neutron-rich nuclei have been 
investigated using theoretical models. Several models of different level of sophistication for determining /3-decay 
half-lives have been proposed and applied over the years. One can mention the more phenomenological treatments 
based on Gross Theory (e.g. lfj) as well as microscopic treatments that employ the proton- neutron quasiparticlc 



random phase approximation (pn-QRPA) [llLll2j or the shell model .13]. The later hybrid version of the RPA models 
developed by Moller and coworkers, combines the pn-QRPA model with the statistical Gross Theory of the first 
forbidden decay (pnQRPA+ffGT) 



141 ] . There are also some models in which the ground state of the parent nucleus is 
described by the Hartree-Fock BCS method, or other density functional method (DF) and which use the continuum 
QRPA (CQRPA) (e.g. Il5|). Recently relativistic pn-QRPA (RQRPA) models have been applied in the treatment of 



neutron- rich nuclei in the N ~ 50, N ~ 82 and Z ~ 28 and 50 regions 



1(5 



17] . Despite continuing improvements the 



predictive power of these conventional models is rather limited for /3-decay half-lives of nuclei. 
The microscopic calculations of allowed weak interaction rates using the pn-QRPA model 
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12, 



18[ and first- 



forbidden /3-decay properties by Homma and collaborators 19] led to a better understanding of the r-process. However 
there was a need to go to domains of high temperature and density scales where the weak interaction rates are of 
decisive importance in studies of the stellar evolution. 

Two different microscopic approaches, namely the shell model 2f]]and the pn-QRPA 21], have so far been succes- 



sively used in the large-scale calculation of stellar weak-interaction mediated rates for r-process applications. Whereas 
the advantage of the shell model clearly lies in it's ability to take into account the detailed structure of the /3-strength 
functions, the shell model is forced to use approximations like Brink's hypothesis (in electron capture direction) and 
back resonances (in /3-decay direction) when extending the calculation to high-lying excited states. These states can 
contribute rather significantly to the total weak rate in high temperature and high density environments. Brink's 
hypothesis states that GT strength distribution on excited states is identical to that from ground state, shifted only 
by the excitation energy of the state. GT back resonances are the states reached by the strong GT transitions in the 
inverse process (electron capture) built on ground and excited states. On the other hand the pn-QRPA model gets rid 
of such approximations and allows a state-by-state evaluation of the weak rates by summing over Boltzmann- weighted, 
microscopically determined GT strengths for all parent excited states. 

In this paper we attempt to calculate the allowed as well as the first-forbidden (FF) /3-decay rates for neutron-rich 
nickel isotopes using the pn-QRPA model for the first time in stellar matter. Our calculations do not employ approx- 
imations (Brink's hypothesis and/or back resonances). The recent measurement of the GT + strength distribution of 



r6 Se 22] supports the argument that, due to nuclear correlations across the N = 40 shell gap. the GT transitions for 
fp-nuclei with proton numbers Z < 40 and neutron numbers N > 40 would not be Pauli blocked. Nevertheless for 
neutron-rich nuclei FF transitions are favored due to the amplification of available phase space. It is in this context 



that we attempt here to calculate both allowed GT and FF transitions for neutron-rich nickel isotopes ( 8 Ni). The 



theory of allowed and first-forbidden /3-decay transitions (along with its leptonic phase-space content and nuclear ma- 
trix elements) is well-established 23W271]. Whereas the allowed /3-decay (having the same parity before and after the 



decay) is relatively simple to calculate, the first-forbidden decay shows a far wider spectrum both in lepton kinematics 



and in nuclear matrix elements 



23 



24, 



27]. 



The pn-QRPA formalism for calculations of the total /3-decay half-lives and FF /3-decay rates for neutron-rich nickel 
isotopes are discussed in section II. Results and comparisons with experimental and other theoretical calculation are 
shown in section III. Section IV outlines the importance of the results in connection with the simulation of the 
advanced stages of stellar evolution and nucleosynthesis and concludes our discussion on consequences of FF /3-decays 
in stellar matter. 



II. FORMALISM 

I - ! 

In the pn-QRPA formalism [28[, proton- neutron residual interactions occur in two different forms, namely as 
particle-hole (ph) and particle-particle (pp) interactions. The particle-hole and particle-particle interactions can be 
given a separable form because P~ transitions are dominated by the particle-hole interaction. These ph and pp forces 
are repulsive and attractive, respectively, when the strength parameters x f° r P n GT force and k for pp GT force 
take positive values. The advantage of using these separable GT forces is that the QRPA matrix equation reduces 
to an algebraic equation of fourth order, which is much easier to solve as compare to full diagonalization of the 



non-Hcrmitian matrix of large dimensionality 
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28 



Essentially we first construct a quasiparticle basis (defined by a Bogoliubov transformation) with a pairing interac- 
tion, and then solve the RPA equation with a schematic separable GT residual interaction. The single particle energies 
were calculated using a deformed Nilsson oscillator potential with a quadratic deformation. The pairing correlation 
was taken into account in the BCS approximation using constant pairing forces. The BCS calculation was performed 
in the deformed Nilsson basis for neutrons and protons separately. The formalism for calculation of allowed /3-decay 
rates in stellar matter using the pn-QRPA model can be seen in detail from Refs. 
the necessary formalism to calculate the FF /3-decay rates. 

Allowed Fermi and GT transitions have a maximum spin change of one unit and no change in parity of the 
wavefunction. In contrast, the most frequent occurrence of a forbidden decay is when the initial and final states have 
opposite parities, and thus the selection rule for allowed decay is violated. The spins of the initial and final states 
can be different at most by AJ= 0, ±1, ±2 for the first forbidden transitions. The isospin selection rule remains the 
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29j . Below we describe briefly 
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same as for allowed decays. For the calculation of the FF /3-decay rates, matrix elements of the separable forces which 
appear in RPA equation are given by 



<W = +2x/pn(M)/ P 'n'(M), (1) 



Cv = -2KfMfv>n>(ti, (2) 

where 

fpn(fi) =< j P mp\t^r[aYi} 2 ^\jnm n > (3) 

is a single-particle FF transition amplitude (the symbols have their normal meaning). Note that /i takes the values 
fi = 0, ±1, and ±2 (for allowed decay rates /i only takes the values and ±1), and the proton and neutron states have 
opposite parities 19[. In the pn-QRPA formalism, proton- neutron residual interactions occur as particle-particle (pp) 
and particlc-hole (ph) interaction. Both pp and ph interaction terms are given a separable form in the model and are 
characterized by two interaction constants (x for the ph force and K for the pp force). For further details we refer to 



19( . Deformations of the nuclei were calculated using 

= 125(Q 2 ) 

1.44(Z)(A) 2 /3' 1 > 

where Z and A are the atomic and mass numbers, respectively and Qi is the electric quadrupole moment taken from 
Ref. [30]. Q-values were taken from the recent mass compilation of Audi et al. 3l|. 



The FF /3-decay rates from the ith state of the parent to the jth state of the daughter nucleus is given by 

Ay — ln2 ^ . (5) 

The phase-space factors are given as integrals over the lepton distribution functions and hence are sensitive 
functions of the temperature and density in stellar interior. The ftij are the comparative half-lives and are related 
to the FF weak interaction matrix elements. 

For the first forbidden transitions the integral can be obtained as 



fij= W» 2 - l(w m -w) 2 [(w m -w) 2 F 1 {Z,w) + (w 2 - l)F 2 (Z,w)](l-G^)dw, (6) 



where w is the total kinetic energy of the electron including its rest mass and w m is the total /3-decay energy 
(w m = m p — mid + Ei — Ej, where m p and Ei are mass and excitation energies of the parent nucleus, and and Ej 
of the daughter nucleus, respectively). G_ are the electron distribution functions. Assuming that the electrons are 
not in a bound state, these are the Fermi-Dirac distribution functions, 

G- = [exp(^^l) + l]-\ (7) 

Here E — [w — 1) is the kinetic energy of the electrons, Ej is the Fermi energy of the electrons, T is the temperature, 
and k is the Boltzmann constant. 

The Fermi functions, F\{±Z, w) and F%(±Z, w) appearing in Eq. (6) are calculated according to the procedure 



adopted by Gove and Martin 



32|. 



The number density of electrons associated with protons and nuclei is pY e N a, where p is the baryon density, Y e is 
the ratio of electron number to the baryon number, and Na is the Avogadro's number. 

^=^<T> a r <G -- G+)!,> ' <8 » 

where p = (w 2 — l) 1 / 2 is the electron or positron momentum, and Eq. (8) has the units of moles cm -3 . This equation 
is used for an iterative calculation of Fermi energies for selected values of pY e and T. 

There is a finite probability of occupation of parent excited states in the stellar environment as a result of the high 
temperature in the interior of massive stars. Weak decay rates then also have a finite contribution from these excited 
states. The occupation probability of a state i is calculated on the assumption of thermal equilibrium, 



= expj-E./kT) 
1 J2 i = 1 exp(-E i /kTY {) 

where Ei is the excitation energy of the state i, respectively. The rate per unit time per nucleus for any weak process 
is then given by 



\=j2 p i^r (io) 



'j 



The summation over all initial and final states are carried out until satisfactory convergence in the rate calculations 
is achieved. 



7 



III. RESULTS AND COMPARISON 

In this section we present the calculation of the total /3-decay half-lives and FF /3-decay rates for neutron-rich Ni 
isotopes. Both allowed GT and FF transitions are considered for calculations of the /3-decay half-lives. The calculated 
half-lives are also compared with the measured ones to estimate the accuracy of the predicted decay rates. 

For the Ni isotopic chains the results are shown in Fig. [1] and compared with experimental values 33H35] and 



with the self-consistent density functional + continuum quasiparticle random phase approximation (DF3 + CQRPA) 

n 

calculation of Ref . 36] . 

The calculation performed in the present work is in better agreement with the experimental data. The inclusion of 
the first-forbidden transitions within the pn-QRPA results in noticeably shorter half-lives and satisfactory agreement 
with experiment is reached when first-forbidden transitions are included. 

The pn-QRPA half-lives are slightly bigger than experimental half-lives 
lives of Ref. 



33 1 for 72 ' 73 Ni whereas the calculated half- 
36j are slightly smaller than the corresponding measured values. For the rest of the cases ( 74_78 Ni) the 



pn-QRPA calculated half-lives are in excellent comparison with the measured half-lives. The calculation shows that 
the inclusion of the FF contribution makes the half-life predictions more reliable for neutron-rich nickel isotopes. 

Tables 1 and 2 display the phase space calculation for allowed and FF transitions, as a function of stellar temperature 
and density, for the neutron-rich nickel isotopes. The phase space is calculated at selected density of 10 2 g/cm 3 , 
10 6 g/cm 3 and 10 10 g/cm 3 (corresponding to low, intermediate and high stellar densities, respectively) and stellar 
temperature Tg = 0.01, 1, 3 and 10 (Tg gives the stellar temperature in units of 10 9 K). Table 1 shows the result for 
72 ~ 75 Ni whereas Table 2 displays the phase space calculation for 76 " 78 Ni. It can be seen from these tables that for 
high densities the phase space increases with increasing stellar temperature. As the stellar core becomes more and 
more dense the phase space decreases. The phase space increases by 1-3 orders of magnitude for FF transitions. This 
enlarged phase space makes the FF /3-decay rates significant for the neutron-rich isotopes as will be discussed below. 

The FF /3-decay rates were calculated for densities in the range 10 5 — 10 11 g/cm 3 and temperature range 0.01 < 
Tg < 30. Figs. [2} [8] show three panels depicting pn-QRPA calculated allowed and FF /3-decay rates at selected 
temperature and density domain for chosen neutron-rich nickel isotopes. The total /3-decay rates are a sum of these 
two contributions. Fig. [2] depicts stellar /3-decay rates for 72 Ni in units of sec -1 . It is pertinent to mention that 
contribution from all excited states are included in the final calculation of all decay rates. It can be seen from this 
figure that for low and intermediate densities (10 2 and 10 6 g/cm 3 ) the /3-decay rates (both allowed and FF) remain 
almost the same as a function of stellar temperature. At low temperatures (Tg < 5) the FF transition is an order 
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of magnitude bigger than the allowed transition primarily because of the enlarged phase space as shown in Table 1 . 
At high stellar temperatures the allowed transitions become dominant as more and more excited states contribute 
significantly to the allowed GT rates. For these high stellar temperatures the allowed GT rates are bigger than the 
FF rates roughly by a factor of five. At high stellar density (10 10 g/cm 3 ) the /3-decay rates tend to go down due to the 
blocking of the available phase space of increasingly degenerate electrons. At high densities and high temperatures the 
contribution to the total /3-decay rates by the excited states is very important. The rates rise with higher temperature 
as more excited states start contributing when the temperature increases. Increasing temperature and density have 
opposite contributions to the calculated rates. Whereas growth of the stellar density suppresses the /3 decay rates due 
to a diminishing phase space available for escaping electrons, increase in temperature weakens the Pauli blocking and 
therefore enhances the contribution of the GT_ transitions from excited states of the parent nucleus. 

Figs. |31 HI and [SJ show the cases of 73 Ni, 74 Ni and 75 Ni, respectively. The results almost follow the same trend as 
shown in Fig. [21 The FF /3-decay rates continue to dominate till still higher temperatures (Tg ~ 20) for the case 
of 74 Ni and 75 Ni (for low and intermediate densities). At high stellar temperatures, Tg ~ 30, the allowed GT rates 
dominate. The effects of increasing temperature and density on the decay rates were discussed earlier. 

Fig. [5] shows a different behavior for the case of 76 Ni. Here the FF rates are bigger than the allowed rates for all 
temperatures at low and intermediate densities. The FF rates are more than an order of magnitude bigger at low 
temperatures. At high temperatures the two rates tend to approach each other with the FF rate still slightly bigger. 
At high densities the FF rates are almost three orders of magnitude bigger at low temperatures. As the temperature 
increases the allowed GT rates tend to increase at a faster rate and at Tg = 30, the allowed rates is almost double 
the FF rate. The reason for this new trend can be traced back partly due to the phase space calculation (see Table 2) 
where the FF phase space is much bigger than the allowed phase space and partly due to calculation of much larger 
values for FF transition amplitudes (nuclear matrix elements). 

Fig. [7] ( 77 Ni) follows the normal trend. One notes that for the case of 77 Ni, for low and intermediate densities, 
the FF rates dominate up till Tg ~ 25. For the case of 78 Ni it can be seen from Fig. [8] that at high densities and 
low temperatures the FF rates are up to two orders of magnitude bigger. As temperature increases the allowed rates 
overtake the FF rates and at Tg = 30 they are almost three times bigger. The reason for this behavior is again 
traced back to phase space calculation for FF transition of 78 Ni (see Table 2). At high densities the phase space 
is around three orders of magnitude bigger for the FF transition and explains the orders of magnitude bigger FF 
/3-decay rates. Increasing the stellar temperature weakly affects the FF /3-decay rates. The reason is that the first 



forbidden transitions are of relatively high energies and basically of the particle-hole nature. Therefore a smearing of 
single-particle strengths due to pairing correlations and temperature influence them considerably less. On the other 
hand the allowed GT rates increase at a faster pace with increasing temperature and hence explains this feature in 
all figures. 

In the /3-decay of a particular isotope, an emitted electron can receive any kinetic energy from zero up to the 
maximum kinetic energy available for the decay. In stellar environments, the electron density can be sufficiently high 
so that the low-energy states of the electron gas are essentially filled. In such an environment, there are fewer states 
available into which a low-energy electron can be emitted, and therefore /3-decay is strongly inhibited if maximum 
kinetic energy available for the decay is small. The /3-decay rate of an isotope depends on the local temperature, the 
local density, the maximum /3-decay energy of the isotope, and the of /3-decay process involved 37[. For FF decays 



(parity changes and AJ= 0, ±1, ±2), when maximum kinetic energy available for the decay is small, the inhibition 
of first-forbidden decays is not as large as the inhibition of allowed, for the same temperature, density, and /3-decay 
energy [36 1. 



IV. CONCLUSIONS 

In this work the total stellar /3 _ -decay half-lives and FF /3-decay rates for a number of neutron-rich nuclei have 
been calculated, for the first time, using the pn-QRPA theory. The temperature and density profiles were relevant to 
the advanced nuclear burning stages of the presupernova and supernova star. The results for total /3-decay half-lives 
are in good agreement with the experimental results 33H35|. 



The microscopic calculation of FF /3-decay rates, presented in this work, could lead to a better understanding of 
the nuclear composition and Y e in the core prior to collapse. These calculated rates are also required in calculations 
of the hydrostatic phases of presupernova evolution through core silicon burning and the onset of core collapse and 
hence are required for an understanding of the explosion mechanism. It has been found that the inclusion of accurate 
weak rates during all stages of nuclear burning, from core oxygen burning onwards, leads to an important difference 
in presupernova iron core structure and its entropy prior to core collapse. 

What impact enhanced FF /3-decay rates may have on core collapse simulation results? The effects of enhanced FF 
/3-decays are two-fold. These enhanced rates on neutron-rich nuclei significantly increase the value of Y e as compared 
to allowed /3-decay rates after oxygen depletion and during silicon shell burning in stellar core. The increment in 
Y e during and after silicon shell burning is primarily due to the increased phase space for the /3-decay electrons and 
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should be more for the reported FF /3-decay rates. Secondly, /3-decays are additional neutrino sources and contribute 
to the cooling of the stellar core and to the reduction in entropy. The lower entropy favors the explosion mechanism. 
The enhanced cooling of the core due to FF /3-decay rates of neutron-rich nuclei will be quite efficient, because the 
average energies of the neutrinos emitted in the FF /3-decays will be bigger than the allowed /3-decay rates and also 
for the competing allowed electron captures. As a consequence, the FF /3-decay rates lead to a significantly lower core 
temperature during and after silicon. 

During the collapse phase allowed /3-decays become unimportant due to the increased electron chemical potential 
which drastically reduces the phase space for the beta electrons. It is well known that a smaller lepton fraction disfavors 
the outward propagation of the post-bounce shock waves, as more overlying iron core has to be photo-dissociated. It 
is noted that the FF /3-decay rates yield higher Y e in comparison to allowed /3-decay rates prior to collapse. This 
might assist the march of post-bounce shock for a successful supernova explosion to occur. These stronger /3-decay 
rates can also assist in a more vigorous URCA process and lead to cooler presupernova cores consisting of lesser 
neutron-rich matter than in presently assumed core-collapse simulations. The allowed and FF /3-decay rates were 
calculated on a fine temperature-density grid, suitable for simulation codes, and may be requested as ASCII files from 
the corresponding author. 

The FF /3-decay rates presented here can be important in determinations of the nucleosynthesis yields of the 
elements and their isotopes, both during hydrostatic burning regimes and the subsequent collapse and explosive 
burning phases. This microscopic calculation of FF /3-decay rates in stellar environment presented here is expected 
to prove quite useful for core-collapse simulator world-wide. 
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TABLE I: Comparison of calculated allowed and first-forbidden (FF) phase space factors for 72 ~ 75 Ni for different selected 
densities and temperatures in stellar matter. The first column gives stellar densities (pY e ) having units of g/cm 3 , where p is 
the baryon density and Y e is the ratio of the electron number to the baryon number. Temperatures (T 9 ) are given in units of 
10 9 K. 





T 9 


Phase space ( 72 Ni) 


Phase space ( Ni) 


Phase space ( 74 Ni) 


Phase space ( 75 Ni) 


Allowed FF 


Allowed FF 


Allowed FF 


Allowed FF 


10 2 


0.01 
1 
3 
10 


1.7E+04 6.5E+04 
4.0E+08 1.0E+10 
4.0E+08 1.0E+10 
4.0E+08 1.0E+10 


9.3E+04 3.9E+05 
2.6E+08 1.4E+10 
2.6E+08 1.4E+10 
2.6E+08 1.4E+10 


8.1E+04 4.4E+05 
1.8E+08 3.2E+09 
1.8E+08 3.2E+09 
1.8E+08 3.2E+09 


4.0E+05 1.7E+06 
5.1E+08 9.9E+08 
5.1E+08 9.9E+08 
5.1E+08 9.9E+08 


10 6 


0.01 
1 

3 
10 


1.7E+04 6.5E+04 
4.0E+08 1.0E+10 
4.0E+08 1.0E+10 
4.0E+08 1.0E+10 


9.2E+04 3.9E+05 
2.6E+08 1.4E+10 
2.6E+08 1.4E+10 
2.6E+08 1.4E+10 


8.0E+04 4.4E+05 
1.8E+08 3.2E+09 
1.8E+08 3.2E+09 
1.8E+08 3.2E+09 


4.0E+05 1.7E+06 
5.1E+08 9.9E+08 
5.1E+08 9.9E+08 
5.1E+08 9.9E+08 


10 1() 


0.01 
1 

3 
10 


0.0E+00 0.0E+00 
5.1E+07 2.3E+09 
5.2E+07 2.4E+09 
6.6E+07 2.8E+09 


0.0E+00 0.0E+00 
1.3E+07 6.6E+09 
1.3E+07 6.6E+09 
2.0E+07 6.9E+09 


0.0E+00 0.0E+00 
9.2E+06 1.1E+09 
9.7E+06 1.1E+09 
1.5E+07 1.2E+09 


0.0E+00 0.0E+00 
4.2E+07 3.2E+08 
4.3E+07 3.2E+08 
5.8E+07 3.5E+08 



TABLE II: Same as Table 1 but for 76 ~ 78 Ni. 





T 9 


Phase space ( 76 Ni) 


Phase space ( 77 Ni) 


Phase space ( 78 Ni) 


Allowed FF 


Allowed FF 


Allowed FF 


10 2 


0.01 

1 

3 
10 


1 qp,_i_f)5 i IK+Ofi 

J. . t/j_j i yjtj i.ii J i uu 

2.3E+08 1.8E+10 
2.3E+08 1.8E+10 
2.3E+08 1.8E+10 


4 w,+05 2 1 E+Ofi 

i . till | \ i ' ) z- . i i t v / v ' 

9.6E+08 1.3E+09 
9.6E+08 1.3E+09 
9.6E+08 1.3E+09 


S 4R+05 2 2E+06 
2.3E+08 4.9E+10 
2.3E+08 4.9E+10 
2.3E+08 4.9E+10 


10 6 


0.01 

1 

3 
10 


2.3E+08 1.8E+10 
2.3E+08 1.8E+10 
2.3E+08 1.8E+10 


9.6E+08 1.3E+09 
9.6E+08 1.3E+09 
9.6E+08 1.3E+09 


2.3E+08 4.9E+10 
2.3E+08 4.9E+10 
2.3E+08 4.9E+10 


10 10 


0.01 

1 

3 
10 


0.0E+00 0.0E+00 
2.2E+07 8.0E+09 
2.3E+07 8.1E+09 
3.0E+07 8.4E+09 


0.0E+00 0.0E+00 
1.3E+08 5.4E+08 
1.4E+08 5.5E+08 
1.7E+08 5.7E+08 


0.0E+00 0.0E+00 
2.5E+07 2.5E+10 
2.5E+07 2.5E+10 
3.3E+07 2.6E+10 
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FIG. 1: Total /3-decay half-lives for Ni isotopes calculated from the pn-QRPA (this work) including only the allowed (GT), 



allowed plus first-forbidden (GT+FF) transitions, in comparison with the DF3+CQRPA 



and experimental data 
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FIG. 4: Same as Fig. H but for 74 Ni 
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FIG. 8: Same as Fig. H but for 78 Ni 



